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ABSTRACT 


A  survey  is  given  of  previous  research  on  the  linear  theory  of  the  lower 
hybrid  drift  instability  (Sec.  I).  A  summary  is  made  of  research  recently  com¬ 
pleted  as  well  as  research  currently  in  progress  on  the  nonlinear  theory  and 


computer  simulations  of  the  lower  hybrid  drift  instability  (Sec.  II). 
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I.  The  Lower  Hybrid  Drift  Instability  -  Linear  Theory 

In  the  past  several  years  there  has  been  considerable  interest  in  the  lower  hybrid  drift  ins¬ 
tability  (Mikhailovskii  and  Tsypin,  1966;  Krall  and  Liewer,  1971),  since  it  may  play  an  impor¬ 
tant  role  in  the  anomalous  transport  processes  of  theta  pinches  and  field-reversed 
configurations.  The  lower  hybrid  drift  instability  may  also  account  for  small  scale  F  region  irre¬ 
gularities  in  the  equatorial  ionosphere  (Huba  et  al.,  1978;  Huba  and  Ossakow,  1979).  The 
existence  of  this  mode  was  first  shown  by  Mikhailovskii  and  Tsypin  (1966)  using  the  cold-fluid 
hydrodynamic  equations.  The  instability  is  caused  by  a  coupling  of  ion  plasma  or  lower  hybrid 
oscillations  to  drift  waves  in  a  nonuniform  plasma. 

Krall  and  Liewer  (1971)  first  treated  this  mode  through  the  Vlasov  equation  with  an 
assumption  kyar«  1  in  a  slab  geometry,  where  a,  is  the  electron  Larmor  radius.  The  usual 
local  approximation  was  used  for  the -profiles  of  the  magnetic  field,  density  and  temperature.  In 
their  model  the  plasma  supports  an  ambipolar  electric  field  E0  in  the  density  gradient  direction. 
It  was  assumed  that  E0  drives  an  electron  current  with  velocity  cEJB0  ,  but  that  the  ions  are 
not  affected  by  the  zeroth  order  electric  or  magnetic  fields.  This  assumption  restricts  their 
analysis  to  frequencies  well  above  the  ion  cyclotron  frequency  and  also  restricts  the  application 
of  their  theory  to  experiments  where  some  feature  prevents  the  electrons  and  ions  from  acquir¬ 
ing  equal  ExB  drifts.  For  example,  in  theta  pinch  experiments,  the  ions  do  not  acquire  a  drift 
because  of  the  lime  scales  involved;  it  is  observed  that  the  time  scale  on  which  the  ambipolar  E 
field  forms  and  excites  the  lower  hybrid  drift  instability  which  then  decays  is  shorter  than  the 
ion  cyclotron  period.  This  same  assumption  has  been  adopted  by  most  of  the  recent  papers. 
They  found  that  the  zero  beta  plasma  is  unstable  to  a  fiutelike  wave,  propagating  across  the 
field,  even  when  T,<  T,  .  Their  condition  for  this  strong  instability  is  that  the  particle  drifts 
and  currents  not  be  too  weak;  vyVf^c,*  ,  where  Vf  is  the  Ex  B  drift  velocity,  v-t  is  the  elec¬ 
tron  diamagnetic  drift  velocity  due  to  the  inhomogeneity  in  the  electron  density,  temperature 
and  magnetic  field,  and  c,  is  the  ion  sound  speed. 


Davidson  and  Gladd  (1975)  dropped  the  assumption  kv at « 1  ,  kept  the  low  beta  limit 
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and  moved  v£  from  the  large-drift-velocity  regime  (  v£»  v„  )  to  the  low-drift-velocity  regime 
(  v£«  V,,  )  where  v„  is  the  ion  thermal  velocity.  Four  additional  important  features  of  the 
lower  hybrid  drift  instability  were  found.  First,  the  wave  has  negative  energy  whenever 
sgn ml ky  —  sgrt  v-t(ky)  .  Hence,  instability  can  exist  even  though  (3/',/8  vy)  v  _„/*  <0  ,  where  F \ 
is  the  ion  velocity  distribution  function.  Second,  the  growth  rate  measured  in  units  of 
,  is  substantial  even  when  v£/v„^l  .  Third,  for  fixed  v£/v„  ,  the  max¬ 
imum  growth  rale  is  an  increasing  function  of  T,  ,  since  the  electron  diamagnetic  current  is 
correspondingly  larger.  Finally,  the  wavelength  that  maximizes  the  growth  rate  generally 
satisfies  kyar=  1  .  For  ,  the  condition  for  y>a>„  requires 

(Vff/2)‘A(m,/m,) 1/4  (Davidson  and  Gladd,  1975;  McBride  and  Hamasaki,  1978),  where 
a,  is  the  ion  Larmor  radius  and  L„  is  the  density  scale  length. 

The  influence  of  finite  plasma  beta  effects  on  the  lower  hybrid  drift  instability,  including 
the  effects  of  transverse  electromagnetic  perturbations  (  SB^O  ),  and  resonant  and  non¬ 
resonant  Vfl0  electron-orbit  modifications,  was  first  investigated  by  Davidson  et  al.  (1976, 
1977).  The  net  effect  of  finite  plasma  beta  is  to  reduce  the  maximum  growth  rate  ym  of  the 
lower  hybrid  drift  instability.  However,  the  details  depend  on  plasma  parameters.  In  the 
regime  where  7",=  T>  and  v£=  v„  ,  finite  beta  electromagnetic  effects  are  destabilizing  for  all  ky 
,  whereas  the  finite  beta  effects  associated  with  resonant  V  Bo  orbit  modifications,  which  lead  to 
a  significant  decrease  in  the  electron  diamagnetic  drift  velocity,  are  stabilizing  for  all  ky  .  The 
two  finite-beta  effects  combine  to  give  a  net  reduction  in  maximum  growth  rate  ym  .  In  the 
cold  electron  (  Te«  T,  ),  low  drift  velocity  regime  (  v£«  v„  ),  the  finite  beta  electromagnetic 
effects  are  destabilizing  for  small  ky  ,  and  stabilizing  for  large  kr  ,  whereas  the  finite  beta  effects 
associated  with  nonresonant  V  Bo  orbit  modifications  are  stabilizing  for  small  ky  and  destabiliz¬ 
ing  for  large  ky.  Combining  these  two  finite  beta  effects  reduces  the  maximum  growth  rate  by 
a  factor  (l+/3,/2)-1/J  relative  to  the  value  obtained  when  /3,-0  .  Finally,  except  in  the  limit  of 
Tf/T,-~Q  ,  a  critical  value  for  the  local  plasma  beta  (  f3cr  )  was  found,  such  that  the  lower 
hybrid  drift  instability  is  completely  stabilized  (  y<0  )  for  /3>/3„  . 
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The  influence  of  magnetic  shear  on  the  lower  hybrid  drift  instability  has  been  siudied 
(Krall,  1977;  Gladd  et  al.,  1977,  Cavioson  et  ai.,  1978;  McBride  and  Krall.  1978),  including 
application  to  microstability  properties  of  a  number  of  experiments.  It  was  found  that 
sufficiently  strong  magnetic  shear  can  completely  stabilize  the  lower  hybrid  drift  instability  for 
kar^\  (Krall,  1977,  Gladd  et  al.,  1977;  Davidson  et  al.,  1978).  However,  computer  simula¬ 
tions  showed  that  shear  reduced  the  growth  rate  and  saturation  level  of  the  instability,  hut  that 
complete  stabilization  of  the  mode  did  not  occur  (Winske,  1978).  A  general  condition  for 
shear  stabilization  of  the  lower  hybrid  drift  modes  was  given  as  L,<  L„(aJ L„+  LJ a,)  for  all 
values  of  v£  and  m/ktv„»\  (Krall,  197"’).  L,  is  the  magnetic  field  scale  length.  Including 
finite  electron  gyroradius  (  ka,^  1  ),  McBride  and  Krall  (1978)  found  that  the  low-density  sys¬ 
tems  (  w pf< ut ct  )  are  much  easier  to  stabilize  with  shear  than  are  higher-density  system  ( 
>• 

All  the  theoretical  work  mentioned  above  have  been  based  on  the  local  approximation  of 
the  gradients  Batchelor  and  Davidson  (1976)  performed  a  nonlocal  stability  analysis  in  which 
the  effects  of  global  equilibrium  properties,  finite  radial  geometry,  and  the  presence  of  a  con¬ 
ducting  wall  were  included  in  their  calculations.  Their  analysis  was  under  assumptions  that  (1) 
the  ion  thermal  velocity  is  much  less  than  the  wave  phase  velocity ,and  (2)  the  characteristic 
instability  wavelength  is  larger  than  a  thermal  electron  Larmor  radius.  In  this  regime,  the  lower 
hybrid  drift  instability  does  not  depend  on  resonant-particle  effects  and  can  be  described  by  a 
macroscopic  model.  They  showed  that  nonlocal  and  local  theories  are  in  agreement  for  the 
fastest  growing  mode.  A  fully  kinetic,  nonlocal,  matrix  dispersion  equation  for  electrostatic 
perturbations  about  a  spatially  nonuniform  cylindrical  plasma  equilibrium  was  derived  by 
Davidson  (1976).  His  analysis  was  carried  out  for  radially  confined  rigid-jotor  equilibria  and 
based  on  the  assumption  of  equilibrium  charge  neutrality.  The  nonlocal  structure  of  the  lower 
hybrid  drift  instability  in  a  reversed  field  configuration  was  investigated  by  Huba  et  at.  (1980) 
by  using  kinetic  theory.  Their  calculation  indudes  electromagnetic  effects  and  V  B  electron 
orbit  modifications,  and  ignores  the  electrons'  ambipolar  £0*5  drift  velocities.  They  found 
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that  the  fundamental  mode  is  well  localized  away  from  the  neutral  line.  Higher  order  modes, 
however,  have  growth  rates  comparable  to  the  fundamental  mode  and  are  much  more  global. 
Recently,  we  have  worked  out  a  nonlocal  analysis  of  the  lower  hybrid  drift  instability  in  the  low 
drift  velocity  regime  by  using  kinetic  theory  and  considering  only  electrostatic  perturbations  in 
slab  geometry.  We  found  that  the  lower  hybrid  drift  instability  is  a  negative  energy  wave  driven 
by  resonant  ions  whenever  sgnut/ky  —  sgn  v-e(ky)  as  predicted  by  the  local  theory.  Further¬ 
more,  it  was  discovered  that  the  lower  hybrid  drift  wave  with  a  finite  kx  value  will  propagate  to 
regions  where  the  electron  drift  velocity  v£  equals  the  wave  phase  velocity  (  o>f  ky  ),  and  be  sta¬ 
bilized  by  these  resonant  electrons.  Our  nonlocal  theory  and  local  theories  are  also  in  agree¬ 
ment  for  the  most  unstable  mode.  Finally,  we  also  found  that  higher  order  modes  have  growth 
rates  comparable  to  the  fundamental  mode  and  are  more  global. 
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II.  Saturation  Mechanisms  and  Transport  of  the  Lower  Hybrid  Drift  Instability 

The  saturation  of  the  lower  hybrid  drift  instability  has  been  studied  by  a  number  of 
authors.  Saturation  mechanisms  include  ion  quasilinear  flattening  (Davidson,  1978),  current 
relaxation  (Davidson  and  Gladd,  1975;  Davidson,  1978;  Davidson  and  Krall,  1977,  Chen  and 
Birdsall,  1980;  Myra  and  Aamodt,  1980),  ion  trapping  (Winske  and  Liewer,  1978;  Chen  and 
Birdsall,  1980),  electron  resonance  broadening  (Huba  and  Papadopoulos,  1978;  Gary  and 
Sanderson,  1979;  Gary,  1980),  electron  trapping  (Drake  and  Lee,  1980)  or  electron  Ex  B  trap¬ 
ping  (Drake  and  Huba,  1981;  Chen  et  al.,  1981)  and  nonlinear  frequency  shift  (Chen  and 
Cohen,  1981;  Ishihara  and  Hirose,  1981).  The  quasilinear  evolution  of  the  lower  hybrid  drift 
instability  was  investigated  by  Davidson  (1978)  for  the  electron  drift  velocities  (  \’d  )  less  than 
the  ion  thermal  velocity.  It  was  shown  that  current  relaxation  (  vd— *0  )  and  plateau  formation  ( 
dFJdvy— 0  )  in  the  ion  velocity  distribution  function  are  generally  competing  processes  for  sta¬ 
bilization.  The  corresponding  saturation  levels  were  given  analytically.  If  the  initial  drift  velo¬ 
city  vd(r~0)  lies  in  the  range 

(45Vn me/8/rt,)l/3v„  ==  vjl  <  vd(t— 0)  <  v„  , 

then  it  may  be  energetically  favorable  for  stabilization  to  occur  through  current  relaxation,  that 
is,  forcing  vrf— 0  .  If 

vrf(/=*0)  <  v„/7  =s  (45\Ztt m,/8 m,) 1/3  v„  , 

then  plateau  formation  will  be  completed  before  vd  relaxes  to  zero.  The  instantaneous  heating 
rates  and  rate  of  momentum  transfer  during  the  current  relaxing  were  estimated  in  Davidson 
and  Gladd  (1975),  Davidson  (1978),  and  Davidson  and  Krall  (1977).  It  was  shown  that  the 
lower  hybrid  drift  instability  can  result  in  substantial  resistivity  and  plasma  heating.  Saturation 
of  the  lower  hybrid  drift  instability  via  current  relaxation  in  the  low  drift  yelocity  regime  was 
observed  by  us  in  our  one  dimensional  particle  hybrid  simulations  when  vd  was  allowed  to  vary 
in  time  self-consistently  during  the  wave  growth  (1980),  and  our  two  dimensional  electrostatic 
particle  simulations  (1981).  Our  two  dimensional  simulations  also  showed  that  plateau  forma¬ 
tion  in  density  rt(x)  occurs  locally  after  saturation  in  the  region  where  electrons  have  the  larg- 
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est  drift  velocities.  Myra  and  Aamodt  (1980)  showed  that  reducing  of  the  relative  eleciron-  ion 
drift  velocity  causes  saturation  of  the  lower  hybrid  drift  instability  in  their  guiding-centers-on- 
axis  model.  However,  Drake  et  al.,  (1981)  showed  that  in  a  finite  beta  plasma,  the  particle 
drifts  and  magnetic  field  are  coupled  and,  consequently,  the  anomalous  dissipation  of  the  parti¬ 
cle  drift  energy  and  magnetic  energy  also  linked.  It  was  found  that  the  magnetic  free  energy 
can  be  substantially  larger  than  the  particle  drift  energy  and  can  effectively  act  as  a  free  energy 
source  to  drive  the  instability.  Therefore,  current  relaxation  will  not  occur  in  the  high  beta 
plasma. 

If  cross-field  drift  velocities  are  kept  constant  in  time,  our  one  dimensional  electrostatic 
particle  hybrid  simulations  (1980)  showed  that  the  lower  hybrid  drift  instability  is  stabilized  by 
ion  trapping.  It  was  found  that  the  saturation  level  predicted  by  quasilinear  theory  (Davidson, 
1978)  is  the  saturation  level  of  the  most  unstable  mode  for  ion  trapping.  In  deriving  the 
saturation  level  for  plateau  formation,  Davidson  began  with  an  energy  conservation  equation 
and  the  only  real  invocation  of  quasilinear  theory  seems  to  be  the  specification  that  saturation 
occurs  when  the  ion  velocity  distribution  function  has  been  "flattened"  around  the  mode  phase 
velocity.  Such  flattening  could  in  principle  be  due  to  a  variety  of  causes  besides  the  usual 
quasilinear  diffusion,  for  example,  trapping.  Furthermore,  he  assumed  that  the  spectrum  is 
sufficiently  peaked  about  kt—kl,  (the  wavenumber  corresponding  to  maximum  growth  for  the 
initial  equilibrium  conditions)  that  l  +  k2/  kH~2  is  a  good  approximation  in  the  integrand  in  his 
calculations.  This  was  equivalent  to  a  single  mode  assumption,  in  which  only  the  most  unstable 
mode  exists.  Finally,  the  trapping  frequency  corresponding  to  the  saturation  level  predicted  by 
quasilinear  theory  is  larger  than  the  growth  rate  and  hence  larger  than  the  bandwidth  Ato=y  as 
well.  Therefore,  ion  trapping  will  be  the  saturation  mechanism  when  the  fastest  growing  mode 
is  dominant,  and  the  saturation  level  predicted  by  Davidson  will  be  the  saturation  level  due  to 
ion  trapping.  Ion  trapping  was  also  observed  by  Winske  and  Liewer  (1978)  in  their  two  dimen¬ 
sional  electromagnetic  particle  simulations  with  larger  than  v„  . 

Electron  resonance  broadening  has  also  been  considered  as  a  possible  saturation  mechan- 
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ism  for  high  beta  plasma  (Huba  and  Papadopoulos,  1978;  Gary  and  Sanderson,  1979;  Gary, 
1980).  Physically,  the  significance  of  this  nonlinearity  is  due  to  the  fact  that  a  substantial 
number  of  electrons  can  attain  a  Vi  drift  velocity  comparable  to  the  phase  velocity  of  the 
wave  in  finite  beta  plasmas.  Over  a  broad  range  of  parameters,  the  saturation  energy  is  less 
than  that  expected  if  stabilization  occurred  through  either  current  relaxation  or  ion  trapping. 
This  stabilization  mechanism  allows  a  steady  state  turbulent  spectrum  to  develop  which  is  cru¬ 
cial  to  several  plasma  phenomena  such  as  magnetic  field  line  reconnection  (Vasyliunas,  1975) 
and  collisionless  shock  waves  (Biskamp,  1973).  However,  this  nonlinear  phenomenon  has  not 
been  observed  in  any  simulation  studies.  This  may  due  to  the  unrealistic  mass  ratio  and  small 
number  of  particles  used  in  simulations. 

In  the  previous  quasilinear  investigation  of  heating  by  the  lower  hybrid  drift  instability, 
only  perpendicularly  propagating  waves  were  considered.  Ions,  which  behave  as  if  they  are 
completely  unmagnetized  (since  <d»(uc,  ),  resonantly  interact  with  the  wave  and  exchange 
momentum  and  energy.  The  electrons,  however,  which  are  lightly  bound  to  the  magnetic  field 
lines  (  ai«to„  )  are  nonresonant  and  therefore  undergo  no  irreversible  energy  or  momentum 
exchange  (neglecting  VB  resonances  which  lead  to  electron  resonance  broadening).  The 
quasilinear  electron  "heating"  previously  calculated  simply  results  from  the  coherent  sloshing  of 
the  electron  velocity  distribution  function  in  the  lower  hybrid  drift  waves  and  is  completely 
reversible.  Drake(1980,  1981)  examined  the  electron  dynamics  in  a  single  large  amplitude,  low 
frequency  (  a>«wCf  )  wave  propagating  perpendicular  to  B  .  At  small  wave  amplitude,  the 
electron  motion  in  the  electric  field  is  simply  given  by  the  usual  Ex  3  and  polarization  drifts 
and  is  accurately  described  by  the  quasilinear  theory.  It  was  found  that  above  a  threshold,  the 
electron  cyclotron  motion  is  strongly  modified  and  the  electron  motion,  becomes  stochastic 
(Drake  and  Lee,  1980).  In  two  dimensional  calculations  (Drake  and  Huba,  1981),  it  was  found 
that  when  electron  Exx  B  drift  velocity  greater  the  wave  phase  velocity,  i.e.,  e<t>!  T,>  (kxL„)~'  , 
this  electron  Ex  B  trapping  can  cause  the  end  of  wave  growth  Then  irreversible  electron  tran¬ 
sport  can  take  place.  An  unstable  lower  hybrid  drift  mode  propagating  along  the  diamagnetic 
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drift  velocity  direction  y  was  found  to  decay  into  off  angle  modes  with  kx~ky  which  saturates 
as  they  £xfl  trap  the  electrons.  The  finite  kx  decay  modes  therefore  act  as  an  intermediary 
which  both  transports  electrons  irreversibly  and  saturates  the  original  instability.  The  off  angle 
modes  with  kx^=ky  have  been  observed  both  in  experiments  (Fahrbach  et  al.,  1979;)  and  in  our 
two  dimensional  particle  simulations  (1981). 

Finally,  saturation  of  the  instability  can  be  due  to  the  frequency  modulation  when  the 
amplitude  of  the  wave  is  not  very  small.  We  (1981)  derived  the  nonlinear  dielectric  response  of 
the  lower  hybrid  drift  instability  in  the  low  drift  velocity  regime  by  solving  a  coupled  Vlasov- 
Poisson  equation  of  a  single  unstable  mode  self-consistently.  The  nonlinear  temporal  evolution 
of  the  mode  was  calculated  analytically.  It  was  found  that  a  finite  perturbation  of  the  ion  orbits 
leads  to  a  nonlinear  frequency  shift  that  reduces  the  mode  frequency  and  has  a  weak  stabilizing 
effect  on  the  lower  hybrid  drift  instability.  The  nonlinear  frequency  shift  does  noi  seem  to  be  a 
potent  saturation  mechanism  in  a  collisionless  plasma,  but  may  be  more  relevant  when  there 
are  ion-ion  collisions,  lshihara  and  Hirose  (1981)  derived  the  nonlinear  dispersion  for  the 
hydrodynamic  lower  hybrid  drift  instability  (  vd»  v„  ),  which  is  characterized  by  the  lime 
dependent  complex  frequency  and  the  renormalized  ion  susceptibility.  They  found  that  when 
the  instability  is  saturated  by  the  frequency  modulation,  the  electron  drift  velocity  slows  down 
by  30°/oC>f  its  initial  value. 
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